Local and Parallel Finite Element Algorithms for Eigenvalue Problems

Abstract Some new local and parallel finite element algorithms are proposed and analyzed in this paper foreigenvalue problems.With these algorithms, the solution of an eigenvalue problem on a fine grid is reduced tothe solution of an eigenvalue problem on a relatively coarse grid together with solutions of some linear algebraicsystems on fine grid by using some local and parallel procedure.A theoretical tool for analyzing these algorithmsis some local error estimate that is also obtained in this paper for finite element approximations of eigenvectorson general shape-regular grids.     

AGE method simulations of a turbulent far‐wake compared to spectral DNS

Turbulent flow simulation methods based on finite differences are attractive for their simplicity, flexibility and efficiency, but not always for accuracy or stability. This paper demonstrates that a good compromise is possible with the advected grid explicit (AGE) method. Starting from the same initial field as a previous spectral DNS, AGE method simulations of a planar turbulent wake were carried out as DNS, and then at three levels of reduced resolution. The latter cases were in a sense large‐eddy simulations (LES), although no specific sub‐grid‐scale model was used. Results for the two DNS methods, including variances and power spectra, were very similar, but the AGE simulation required much less computational effort. Small‐scale information was lost in the reduced resolution runs, but large‐scale mean and instantaneous properties were reproduced quite well, with further large reductions in computational effort. Quality of results becomes more sensitive to the value chosen for one of the AGE method parameters as resolution is reduced, from which it is inferred that the numerical stability procedure controlled by the parameter is acting in part as a sub‐grid‐scale model. Copyright © 2002 John Wiley & Sons, Ltd.     

Multilevel preconditioners for a quadrature Galerkin solution of a biharmonic problem

Efficient multilevel preconditioners are developed and analyzed for the quadrature finite element Galerkin approximation of the biharmonic Dirichlet problem. The quadrature scheme is formulated using the Bogner–Fox–Schmit rectangular element and the product two‐point Gaussian quadrature. The proposed additive and multiplicative preconditioners are uniformly spectrally equivalent to the operator of the quadrature scheme. The preconditioners are implemented by optimal algorithms, and they are used to accelerate convergence of the preconditioned conjugate gradient method. Numerical results are presented demonstrating efficiency of the preconditioners. © 2005 Wiley Periodicals, Inc. Numer Methods Partial Differential Eq 2006     

Finite Element Solution of Conical Diffraction Problems

This paper is devoted to the numerical study of diffraction by periodic structures of plane waves under oblique incidence. For this situation Maxwell's equations can be reduced to a system of two Helmholtz equations in R ² coupled via quasiperiodic transmission conditions on the piecewise smooth interfaces between different materials. The numerical analysis is based on a strongly elliptic variational formulation of the differential problem in a bounded periodic cell involving nonlocal boundary operators. We obtain existence and uniqueness results for discrete solutions and provide the corresponding error analysis.     

A finite model property for RMImin

It is proved that the variety of relevant disjunction lattices has the finite embeddability property. It follows that Avron's relevance logic RMI _min has a strong form of the finite model property, so it has a solvable deducibility problem. This strengthens Avron's result that RMI _min is decidable. (© 2006 WILEY‐VCH Verlag GmbH & Co. KGaA, Weinheim)     

Indirect Methods For Determination of The Protective Effects of Coating Films on The Surface of Crystals

The extents of the protective effects of coating films on the surface of crystals were determined. Three different samples were made with different quantities of coating fluid (Sepifilm LP 010 in 10% aqueous solution). Since the atomizing rate was constant, the coating time increased in parallel with the volume of coating fluid applied. The direct measurement of film thickness and smoothness is very difficult, and therefore indirect methods were used. Dimenhydrinate was chosen as model drug; this is a heat-sensitive antihistamine with a low melting point. This temperature can be reached during the tableting process. The behaviour of samples on exposure to heat was examined by differential scanning calorimetry. The water uptakes of the samples were determined with an Enslin apparatus. Plasticity was studied with an instrumented tablet machine. These indirect methods (thermal conductivity, water uptake and plasticity measurements) revealed connections between the results of the various experiments. An overlong coating time decreased the protective effect of the coating film. This revised version was published online in August 2006 with corrections to the Cover Date.     

Sub(L)中的五边形格特征及其数量不等式

五边形结构在刻画格的特征方面具有十分重要作用，应用格论及组合数学的方法讨论了格L与其子格格Sub(L)中所含五边形格之间的数量关系，给出了有限格L的子格格中三个元生成五边形格的充要条件，同时给出了Sub(L)所含不同五边形格数量的一个下界。     

Localization effects and measure source terms in numerical schemes for balance laws

This paper investigates the behavior of numerical schemes for nonlinear conservation laws with source terms. We concentrate on two significant examples: relaxation approximations and genuinely nonhomogeneous scalar laws. The main tool in our analysis is the extensive use of weak limits and nonconservative products which allow us to describe accurately the operations achieved in practice when using Riemann-based numerical schemes. Some illustrative and relevant computational results are provided.

     

高折射率椭圆芯布拉格光纤的偏振特性

应用带完全匹配层边界条件的全矢量有限元方法,分析了高折射率椭圆芯布拉格光纤的偏振特性,详细讨论了光纤结构参量对模式双折射度以及群速度走离的影响,结果发现:高折射率椭圆芯布拉格光纤的模式双折射度可达10-2量级,比传统保偏光纤至少高出一个量级,并且表现出不同于传统保偏光纤的群速度走离特性;在保持较高双折射度的同时,通过合理的结构设计,可在一定的波长处灵活地获得较大的群速度走离或零走离特性,具有重要的潜在应用价值;最后,简要分析了低折射率区域的折射率变化对偏振特性的影响;研究结果有助于设计高性能的微结构保偏光纤。     

Semi‐analytical integration of the 8‐node plane element stiffness matrix using symbolic computation

The semi‐analytical integration of an 8‐node plane strain finite element stiffness matrix is presented in this work. The element is assumed to be super‐parametric, having straight sides. Before carrying out the integration, the integral expressions are classified into several groups, thus avoiding duplication of calculations. Symbolic manipulation and integration is used to obtain the basic formulae to evaluate the stiffness matrix. Then, the resulting expressions are postprocessed, optimized, and simplified in order to reduce the computation time. Maple symbolic‐manipulation software was used to generate the closed expressions and to develop the corresponding Fortran code. Comparisons between semi‐analytical integration and numerical integration were made. It was demonstrated that semi‐analytical integration required less CPU time than conventional numerical integration (using Gaussian‐Legendre quadrature) to obtain the stiffness matrix. © 2005 Wiley Periodicals, Inc. Numer Methods Partial Differential Eq, 2006